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Abstract 

We show that a double seesaw framework for neutrino masses with fi — r exchange symmetry 
can lead to one of the righthanded seesaw partners of the light neutrinos being massless. This 
can play the role of a light sterile neutrino, giving a 3 + 1 model that explains the LSND results. 
We get a very economical scheme, which makes it possible to predict the full 4x4 neutrino mass 
matrix if CP is conserved. Once CP violation is included, effect of the LSND mass range sterile 
neutrino is to eliminate the lower bound on neutrinoless double beta decay rate which exists for the 
three neutrino case with inverted mass hierarchy. The same strategy can also be used to generate 
a natural 3 + 2 model for LSND, which is also equally predictive for the CP conserving case in the 
limit of exact u — r symmetry. 
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I. INTRODUCTION 



Our knowledge of neutrinos has undergone a major revolution in the last seven years. 
There is now convincing evidence for nonzero neutrino masses and mixings. Most of the 
data seem to be consistent with a picture of three active neutrinos (is e , u^, v T ) having mass 
and mixing among themselves. While two of the three mixings that characterize the three 
neutrino system i.e. #12 corresponding to solar neutrino oscillation and #23 corresponding 
to that for atmospheric neutrinos are fairly well determined, the third mixing angle #13 has 
only an upper limit [1] and is the subject of many future experiments 

as. 

In anticipation of the results from the planned experiments, there has been theoretical 
speculations about possible meaning of the existence of a small or "large" ^13 j3] - An inter- 
esting possibility which has been inspired by the experimental observation of near maximal 
atmospheric mixing angle (#23) and small upper limit on (9 13 is that there may be an ap- 
proximate 11 — r interchange symmetry in the neutrino sector |5j]. Even though large mass 
difference between the muon and the tau would appear to go against it, in supersymmetric 
theories where the muons and the neutrinos get mass from different Higgs fields, we may 
have a situation where an exact fj, — r symmetry prevails at high scale but at lower scales 
it may only survive in the H u sector but not in the Hd sector. Examples of this type exist 
in literature 0, Q|. The next generation of #13 searches can throw light on this symmetry 8] 
making it an important idea to pursue. Approximate u — r symmetry has also interesting 
implications for the origin of matter via leptogenesis{2|. In this paper, we discuss another 
application of the idea of /i — r symmetry for neutrinos. 

While most of the present data have been discussed in the context of a three neutrino 



picture, the results of the LSND experiment [lo|. if confirmed by the Mini-BooNe[ll| experi- 
ment currently in progress would require that there be one 3, 3 or two sterile neutrinos 
that mix with the active ones. The sterile neutrinos must have a mass in the range of an eV, 
thereby posing a major challenge for theories since being standard model singlets the sterile 
neutrinos could in principle have a mass of order of the Planck scale. There exist a num- 
ber of interesting proposals in literature to explain the lightness of the sterile neutrino 
and in this paper we show that this can happen in the context of the conventional see- 
saw mechanism[16] which is anyway needed to understand the small masses of the known 
neutrinos. 
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Furthermore, with four neutrinos mixing with each other, it becomes difficult to make any 
predictions for the various mass matrix elements. Our model has the additional advantage 
that in the CP conserving case, the symmetries of the theory allow us to write down all 
the elements of the 4x4 mass matrix. For some examples where predictions for the four 



neutrino case have been attempted, see 
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The seesaw mechanism introduces three right handed neutrinos to the standard model. 
Since the right handed neutrinos are sterile with respect to known weak interactions, it 
has often been speculated as to whether one of them can somehow be ultralight and play 
the role of the sterile neutrino while at the same time maintaining the conventional seesaw 
mechanism for the active neutrinos. Clearly, the only way it can happen is if one can "turn 
off" the expected large mass of one of them without effect of the seesaw mechanism. The 
first example of a model where such a situation arises was proposed in Ref-fl^ where the 
neutrino sector was assumed to obey an (L e — — L T ) <g> S2R symmetry. In this case, one of 
the three right handed neutrinos becomes massless instead of being supermassive as in the 
usual seesaw case while the other two become superheavy. The latter become responsible 
for the smallness of the masses of two of the active neutrinos via seesaw. The third one does 
not need a superheavy seesaw partner since it also remain exactly massless as a consequence 
of the symmetry. The effect of the symmetry is therefore essentially to turn the 3x3 seesaw 
into a 2 x 2 one. Then turning on small deviations from this symmetry one can generate 
tiny mass for the v a as well as small mixings with v e and v u . Using this idea, one could write 

n n 

models for both 3+1 19] as well as 2+2 [20J scenarios. 

In this paper, we point out that in the context of a generalized seesaw model, a somewhat 
simpler symmetry <8> S^r, leads to a predictive 3 + 1 picture for the LSND experiment 
which is consistent with all observations. As noted earlier, there are additional motivation 
to consider such symmetries. 

We assume that the usual seesaw mechanism or more precisely the right handed neutrino 
mass matrix originates from an effective field theory whose ultraviolet completion contains 



the double seesaw form of the neutrino mass matrix [2l( operating above the seesaw scale. 
The Sfj, T) L <£> S^r is assumed to operate at this high scale. We then assume that there is a 
low scale symmetry breaking which gives mass to the sterile neutrino as well as its mixings 
with the active neutrinos. We find that in the absence of CP violation, as long as S^t.l 
symmetry is exact, we can predict the full 4x4 neutrino mass matrix. The active neutrinos 
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have an inverted hierarchy. We further point out that once we include CP violation, even 
though active neutrino hierarchy is inverted, there is no lower bound on the effective mass 
probed by the (5(3q v decay. 

We also show that the same strategy can be extended to generate a natural 3 + 2 model 
for LSND where one of the sterile neutrino is the right handed seesaw partner and another 
is a standard model singlet employed to implement the double seesaw scheme. In the same 
way as in the 3 + 1 case, for the CP conserving case, here also can we predict all elements 
of the 5x5 neutrino mass matrix using existing fits to oscillation data. 



II. THE BASIC IDEA AND THE S^L ® S^R MODEL 

The basic idea for getting the 3+1 model from the type I seesaw was outlined in Ref . jl9^ . 
The seesaw formula for light neutrino mass matrix is given by: 

M v = -MlM^Mo (1) 

Suppose we have a symmetry S in the theory that leads to DetMn = with only one of 
the eigenvalues zero. If the same symmetry also guarantees that DetM^ = 0, then one can 
use the seesaw formula. The way to proceed is to "take out" the zero mass eigen states from 
both Mr and Md and then use the seesaw formula for the 2x2 system. It is then clear 
that the spectrum will consist of two light Majorana neutrinos (predominantly left-handed), 
two heavy right handed Majorana neutrinos; one massless right handed neutrino, which will 
play the role of the sterile neutrino of the 3+1 model and a massless left-handed neutrino. 
Breaking the symmetry S very weakly or by loop effects can lead to a nonzero mass for 
the sterile neutrino as well as its mixings with the v e ^ so that one has an explanation of 
the LSND result via v e — > v s — > process. As already noted in Q], S was chosen to be 
(L e — L M — L T ) ® S 2 r- Here we take a different path. 

We consider a model based on the gauge group SU (2) l x U(1)i 3R xU(1) b-l ■ In addition 
to the standard model matter we have SU(2)l singlets x C ; Nf and Si (i = e, /i, r) with 
quantum numbers given in table 1 (we do not bring in the quark sector, whose gauge 
quantum numbers are obvious) . Here we have used the supersymmetric chiral field notation 
and dropped the chirality (L,R) indices. The superpotential involving these fields can be 
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written as follows: 



W 



KLH U N C + fN c X c S + \<j)SS 



(2) 



The couplings h v , f, A are all 3x3 matrices in flavor space whose properties are constrained 
by the symmetry S^l <E> S^r as discussed below. 



field 


SU(2) L 




U(1) B -l 


Li = {y, e)i 


2 





-1 


N? 


1 


i 

~2 


1 


Si 


1 








H u 


2 


1 

2 





H d 


2 


1 

2 





x c 


1 


1 
2 


-1 





1 









Under S^l symmetry, the left handed leptonic doublets (L^,L T ) transform into each 
other; similarly under S^r, (N°,N°) transform into each other. Invariance under this 
symmetry then restricts the matrices M D = h u v u and M = f < x c > as follows: 

/ mn rn 12 m 12 \ 



M, 



D 



m 2 i m 2 2 m 22 
\ m 21 m 22 m 22 ) 



(3) 



and 



/Mn M 12 M 12 \ 
M = M 21 M 22 M 23 (4) 
V M 21 M 22 M 23 / 

The matrix fj, on the other hand has an arbitrary symmetric form. Our considerations 
are independent of whether the mass matrix elements are real or complex. We keep them 
complex to be quite general. 

The superpotential of Eq. (2) leads to the following double seesaw form for the neutrino 
mass matrix at the scale \x = A < <fi >^> M ^> v w k- 

( m D 



M 



uNS — 



m T D 



V o 

5 



Ml 



\ 

M R 

H J 



(5) 



The zeros in this mass matrix can be guaranteed by an extra U(l)x symmetry under which 
the fields iV c , L, S, 4> and \ c have charges 1, —1, 2, —4 and —3 respectively. Note that we do 
not include a (ft superfield in the theory. 

In order to obtain the light neutrino masses and mixings, we proceed in two steps: first, 
we decouple the S-fermions whose masses are above the scale M and write down the effective 
theory below that scale. The effective u, N mass matrix below this mass scale can be written 
as: 

/ m D \ 
\m T D M^M 1 J 

From Eq. (4), we see that determinant of the (N e , N^, N T ) mass matrix Mr = M/i _1 M T 
vanishes and the zero mass eigenstate is given by 775 (A^ — N T ). We "take out" this state 
and the remaining right handed neutrino mass matrix is 2 x 2 involving only the states 
(N e , 771 (A^ + N T )). Similarly, we can take out the state 4g (u^ — u T ) from the Dirac mass 
matrix nip again leaving a 2 x 2 mass matrix in the basis as 

1 f N e \ 

As already noted the seesaw matrix is now a two generation matrix giving two light Majorana 
neutrino states. The other active light neutrino state is the massless state ^7= (z/ M — u T ). The 
light sterile neutrino state is the massless state v s = (N^ — N T ). This therefore provides 
the starting ingredient of the 3+1 sterile neutrino model. 

III. SYMMETRY BREAKING, MASSES AND MIXINGS OF THE STERILE 
NEUTRINO 

Having explained why the sterile neutrino is light, we now proceed to make it a realistic 
model for LSND by putting a mass term of order eV for v s and smaller off diagonal mixings 
for v s with the active neutrinos. These small scales could arise from a scalar field having a 
keV scale vev as in Ref. 



which can also help us to avoid the WMAP bound 



So in the basis (v e , u+, i/ s ), we have a mass matrix of the form 

M v , v „c& = m af3 (8) 

There are two cases that one can consider: (i) only S^r is broken so that the massless state 
z/_ remains decoupled and we have a 3 x 3 mass matrix with involving (u e , u + , u s ) and (ii) 



when all four states are mixed and both the interchange symmetries are broken. The second 
case is similar to a recent analysis in Ref . . 



A. Case (i): unbroken S^^l 

The light neutrino mass matrix in this case is a 4 x 4 matrix which in weak (flavor) basis 
can be written as 

(m n mia m 12 f x \ 
mi 2 m 22 m 22 h 
m 12 m 22 m 22 f 2 
V fi h h fJ 
Equality of certain entries in this matrix is a reflection of the unbroken symmetry. This 
matrix in general has six real parameters and three complex phases. To diagonalize this 
matrix, we first note that without mixing terms fi, f 2 , the diagonalizing orthogonal matrix 
is given by 



(9) 



( c s 0\ 

— — — n 

V2 V2 V2 U 

— — — n 

V2 V2 V? 



V o o 01/ 

where s = sin 9, c = 0,and 9 is solar mixing angle and eignevalues 

1 



(o) mn + 2m 22 , , 
mr = -\/ 8m i2 + (2m 22 - mn) 2 



( o) mn + 2m 22 
m 2 = 



+ ^\/ 8m 12 + ( 2m 22 - W-ll) 2 



771 



(0) 



(0) 





/4 



(10) 



(11) 



Note that the vanishing of the mass of the third eigenstate is due to unbroken S^x symmetry. 
This case has inverted hierarchy for active neutrino masses. 

In the limit where the sterile neutrino decouples, fj, — r symmetry would have dictated 
that we have four parameters describing the 3x3 neutrino mass matrix. However, in our 
case, the effective theory has a higher symmetry M v = M u S IJiTi l = S^^My. As a 
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result, since there are three experimental inputs i.e. Am 2 &1 Am 2 atm and sin 2 # , the three 
neutrino mass matrix in our case is completely fixed by data. Using the best fit values 



Am 2 



2.2 x 10~ 3 eV 2 , Ami = 7.9 x 10~ 5 eU 2 and sin 2 = 0.30 



active neutrinos is given by 



24j,the mass matrix of 



(12) 



/ 0.0471547 0.000270318 0.000270318\ 
M v = 0.000270318 0.0237444 0.0237444 (eV). 
V 0.000270318 0.0237444 0.0237444 J 

We have ignored CP violation in this discussion. From Eq.(12), we can read off the value 
for the effective neutrino mass < m ee > probed by (3[3ov experiment. It is nothing but the ee 
entry in Eq.(12) and we find that the central value of < m ee >~ 47 meV. In this decoupled 
case, of course, we cannot accomodate the LSND results since the decoupling essentially 
amounts to very small mixings. 

Including the effect of the small terms fx, f 2 <C f± so that LSND mixings are accomodated 
in the picture as originally intended, leads to the following form for the generalized 4x4 
PMNS matrix. 



U 



( 


c 


s 









— s 
V2 


c 

V2 


i 

V2 

-1 






— s 
V2 


c 
V2 


$2 


V 


— o\c + y/28 2 s 


—5is — V25 2 c 





1/ 



(13) 



where 



fx f /- 



x J 1 x 

0i - -r;°2 
H 



2 

h 



(14) 



f 2 f 2 

The mass eigenvalues now get shifted to the values m\ — -h and m^—f^. LSND results require 

f f 2 
that ^ ~ 0.2. Naturalness would require that ~ 0.01 implying that fx >2 — 0.05 eV. A 

detailed numerical analysis of this which fits the central values of all observed parameters is 

given below. 

Unlike the case when the sterile neutrino is decoupled, here we have six parameters 
describing the four neutrino mass matrix in the absence of CP violation. All the elements 
- *. _ m a trix can th e refore b e det e rmined Using - , a nd ^ fr „ m Q - 



8 



well as /S.m 2 LSNDl \U e ±\ and |?7 c/ ,||l4j. We have used /S.m 2 LSND = 0.92eV 2 , |£7" e4 | = 0.136 and 
\Ufj,4:\ = 0.205 and find (ignoring CP phases) that the four neutrino mass matrix is given by 

/ 0.06214 0.0240057 0.0240057 0.117 \ 



(eV) 



(15) 



0.0240057 0.0612706 0.0612706 0.18441 
0.0240057 0.0612706 0.0612706 0.18441 
V 0.117 0.18441 0.18441 0.955 / 
We see that in the absence of CP violation, the value of < m ef _>'^. 62 meV. This is accessible 



to the proposed searches for neutrinoless double beta decay 



25] 



B. Case (ii) 

We can extend U to the more general case where the S^^l symmetry is also broken. In 
this case, the neutrino mass matrix has all entries arbitrary except that / 4 is the largest 
entry in it. Since is broken weakly, the U matrix can be written as 



U 



— s 

— 8 

v/2 



V2 
c 

V2 

~ 6lS — 



013 

1 

V2 
-1 
V2 

V2 



U 



(16) 



613 in this case will be proportional to the small S^ t ,l breaking. However, in this case, even 
in the absence of CP violation, one cannot predict the neutrino matrix fully. 

IV. PHENOMENOLOGICAL IMPLICATIONS 



An interesting implication of the sterile neutrino is its impact on the neutrinoless double 
beta decay. As we saw in the previous section, if CP phases are completely ignored, we get 
a definite prediction for the effective double beta decay mass < m ee > which is the ee entry 
of the neutrino mass matrix. Once CP violation is included, the situation changes and we 
have the effective electron neutrino mass given by 



< m ee >~ \c z mi + s 2 e 2ia m 2 + |C4 4 Pe^m 4 



Am 2 atm (c 2 + s 2 e 2ta ) + \U eA \'e 



2 a 2if3J A 2 

iAm LSND 1 



(17) 

where Am^ m = m| and Am 2 LSND ~ ml are the atmospheric and the LSND mass squared 
differences . Note that which < m ee > can be large (if there is no accidental cancellation as 



was the case in Eq. (15)) with 



< m RR > 



AmL + \ U e4 \ 2 



<^ m LSND — 



62 meV 



(18) 



In figure 1 we present a scatter plot for < m ee > as function of |£/ e4 | 2 ~ 8f. One can see that 
there is a parameter range where the rate of (/3/3)o„ process almost vanishes, which implies 
that there is CP violation due to sterile neutrino sector. 

100r 




FIG. 1: Scatter plot of \m e ff\ =< m ee > in the case of inverted hierarchy in the presence of a 
sterile neutrino with mass ~ 1 eV as a function of 5f. Note that the minimum value of \m e ff\ 
depends on the best fit value for the LSND mixing parameter in a 3+1 scheme. The scatter of 
points is due to the variation of the phases a and (3 between to ir. 



Thus the well known lower bound on neutrinoless double beta decay in the case of in- 
verted hierarchy [26] disappears once the effect of sterile neutrino is included. In other words, 
if Mini-BooNe confirms LSND results, a null result in the next round of double beta decay 
experiments will not necessarily rule out the inverted hierarchy scenario for neutrinos nor 
the Majorana nature of the neutrino. For example, in the absence of the sterile neutrino, if 
there was no signal for f3f3o v decay at the level of 30 meV and if long base line experiments 
confirmed that neutrinos have inverted hierarchy, one would have concluded that the neu- 
trino could be a Dirac fermion. The presence of the sterile neutrino would clearly change 
this conclusion. This is an important effect of the presence of the LSND sterile neutrino. 

We also note that the neutrino mass inferred from the study of the shape of the electron 
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energy spectrum near the end point of radioactive nuclear decay is 



5>. 



ei I 



AmL + \U e4 \ 2 \/Am 2 LSND ~ 62 meV 



(19) 



This is below the present expectations in the Katrin experiment 
future studies in this area. 



271 ] but is of interest for 



V. EXTENSION TO 3 + 2 CASE 



My 



NS 



(20) 



In this section, we discuss how our strategy can be extended to obtain a 3 + 2 model 
for LSND. There are several ways to achieve this. The essential new step is to extend the 
discrete symmetry to S^t.l x S^r x <S^t,s an d introduce an additional singlet fermion. This 
implies that in the symmetry limit only the combination S 2 + S3 field has mass and there 
are two massless sterile neutrinos: N° — N° and — S T . We can now give mass to them 
as well as mixings by including the Majorana mass terms of type i^S a Sb where \i is inthe 
eV range, that break S^r X S^s part of the discrete group while keeping S^l unbroken. 
The form of the mass matrix in this case is given by 

/ m D 
ml M RS 
M T RS 
V M ST 

This reduces to the form in Eq. (5) as /i becomes large giving two light active neutrinos via 
the double seesaw as discussed previously. 

With 2 sterile neutrinos with \i — t symmetrythe general mass matrix can be written as 

/ mu ^12 m 12 fi 9i \ 
mi2 m 2 2 m 2 2 h 92 
m 12 m 2 2 m 2 2 h 92 
h h h 9s 
V 9i 92 92 93 9a J 
The parameter g 3 can be set to zero by a rotation of the two sterile neutrinos leaving a total 
of nine parameters. Compared with 3+1 scheme, we have three extra parameters, gi,g 2 ,g4 : - 



\ 


M ST 
» ) 



M,, 



(21) 
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From the experiment data fit, we have Am^, \U e i\, \U^\, Am^, \U e s\, |L^s|[14], which has 
three more inputs, we can determine the full mass matrix (in the absence of any CP phases). 
Now the 5x5 PMNS matrix U is given by 



U 



c 


s 





Si 


61 \ 


—s 
V2 


c 
V2 


1 

V2 


5 2 




— s 
V2 


c 


-1 
v/2 


82 


£2 


+ V25 2 s 


—5is — V2S2C 





1 





+ V2e 2 s 


— eis - y/2e 2 c 








1 J 



(22) 



All the elements of the mass matrix can be determined using Ami,, Am 2 atm and sin 2 6 ^ as 



well as Am^, \ U e 4\, \ U^\, Am^, \ U e5 \, |?7 r t;||l4j. Corresponding to the two solutions in 
we get two solutions for the mass matrix. 

With the best fit data (i) Am 2 , = 0.92eV 2 ,\U e ±\ = 0.121, {U^ = 0.204, Am 2 51 
22eV 2 , \U e5 \ = 0.036, \U^\ = 0.224|lJ. 

/ 0.0672934 0.0618002 0.0618002 0.110382 0.167045 \ 

0.0618002 0.299067 0.299067 0.186183 1.04006 

M v = 0.0618002 0.299067 0.299067 0.186183 1.04006 

0.186183 0.964982 0.00456411 

1.04006 0.00456411 4.69549 / 



da 



0.0618002 0.299067 
0.110382 0.186183 
1.04006 



V 0.167045 

e 

0.89eF 2 , \U e5 \ = 0.125, \U^\ = 0.160 



(eV) (23) 



With the best fit data (ii) Am 2 41 



0A6eV 2 , \U ei \ = 0.090, \U. 



/ill 



0.226, Am 2 51 



(eV) (24) 



/ 0.06742 0.0329899 0.0329899 0.0568206 0.11209 \ 

0.0329899 0.0826493 0.0826493 0.14289 0.143498 

-U>= 0.0329899 0.0826493 0.0826493 0.14289 0.143498 

0.0568206 0.14289 0.14289 0.685108 0.00398793 

V 0.11209 0.143498 0.143498 0.00398793 0.947753 / 

In both the cases, there is a prediction for m ee ~ 67 meV, which is very similar to the 
3 + 1 case without CP violation. Once CP violating phases are included, there is an effect 
on m ee . The effects for both cases are again similar to the 3 + 1 case. 
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VI. SUMMARY 



In summary, we have presented an example where using a simple discrete symmetry, one 
gets one of the seesaw right handed neutrinos to play the role of a light sterile neutrino 
needed to understand the LSND results. The discrete symmetries used are experimentally 
motivated by the near maximal value of the atmospheric mixing angle. The model not 
only can accomodate the LSND results, but we find that due to the symmetries used to 
generate the light sterile neutrino, all elements of the neutrino mass matrix are determined 
by present data if CP phases are ignored. Once CP violation is included, the lower bound 
on the effective mass for neutrinoless double beta decay known for the three neutrino case 
with inverted hierarchy is eliminated. Thus confirmation of LSND results by Mini-BooNe 
experiment will not only revolutionize our thinking about the nature of new physics beyond 
the standard model (by requiring the concepts such as mirror neutrinos or leptonic symme- 
tries or something else) but it will also affect our conclusions about whether the neutrino is 
a Majorana or Dirac fermion. We further show that our strategy and analysis can be very 
simply extended to the 3 + 2 scenario for LSND, with similar results. 

We would like to thank J.Schechter for useful discussions. This work is supported by the 
National Science Foundation grant no. Phy-0354401. 
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